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H,N» Krltfkos 

Moore School of Electrical Engineering 
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ABSTRACT 


It Is shown that In band-limited (or visible angle limited) 
extrapolation the larger effective aperture L' that can be realized 
from a finite aperture L by over sampling is a function of the accuracy 


of measurements 

1*1 <i . 


. It is shown that for sampling in the interval 
b > 1 the signal must be known within an error e^ given 


by 


2 

®N 


1 

4(2kL») 3 


_e_ JL 
8b L' 


2kL’ 


where 

L is the physical aperture 
L* is the extrapolated aperture 
k ■ 2ttA > A 1b the wavelength 


1. 


BACKGROUND 


(I 


The extrapolation problem finds Its origin in the original in- 

1 2 

vest lgat Ion reported by D. Steplan and H. Pollack 9 on the Spheroidal 
functions and the results on the Uncertainty Principle reported by 
H. Landau and H. Pollack. Recently a resurge of interest in this 
area has resulted in large numbers of published papers two of which 

t i« d 

written by A. Jain and R, Schafer give an excellent review of the 
scope of the present research. K. Abend 7 also has given a review of 
the limitations of extrapolation methods. 

One of the recurring problems occurs in any attempt to Implement 
a useful algorithm of extrapolation is the numerical instability. Since 
the problems are nearly ill posed small errors are capable of generating 
totally erroneous results. In this communication an attempt will be 
made to estimate the bounds for the acceptable errors and then place a 
celling of the measurement accuracy and computational accuracy needed 
for the extrapolation. 


2. PROBLEM STATEMENT 


Consider the function F(x) of finite norm known in the interval 
-L/2 < x < L/2, The Fourier transform F(u) is band limited in the in- 
terval I u I < 1 • One has the well known relation 

’ o 



h(x - x’> g(x*> dx 1 


F(x) 


( 1 ) 
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«» * ■ 




The problem of interest here is ee follows: 

Given the function F(x) in the internal -L < x < L find the source 
function g(x') end then extrapolate into the region |x| > L. In its 
corresponding discrete fora Steplan has shown that although in general 
a non unique answer exist the requirement that the solution has a minimum 
norm ensures the uniqueness of the solution. 

Despite the simplicity of the formulations numerical solutions are 
very difficult or almost impossible to obtain* The gCx*) Is wildly ■vary*" 
lng function as can be Inferred from the following argument. 

In the space of spheroidal orthogonal functions ^ n (x) one observes 
the following; 


F(x) - ld n i|> n (x) , ||F<x)|| - za n < « (2) 


and 


h(x-x') ■ I ij> n (x) ^<x’) * ||^ n (x)||-l 

a a 2 

g(x’> - Z ^ n Cx’) , ||g(x)|| - E -j- 2 - 


(3) 

C4) 


where are the eigenvalues of the spheroidal functions. The problem 
here Is that for any functions F(x) with bounded norm ||f|| < « the 
source function g(x) may not be bounded. 

The Teaaon is that when one attempts to extrapolate beyond the 
degrees of freedom dictated by the bandwidth and time interval Slepian 
and later H. Landau 6 has shown that the asymptotic form of the eigen- 
values is given by 
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A R Cc) + g- 

R (1 + « 1 


C5) 


when the Index n Is given by 


n . »W «<?) e x b_lgg.,e 
" 2ir „2 

IT 


< 6 ) 


where m(a Q ) it the measure of the u internal and m(T) is the measure of 
a finite set such that [c m(s)] - 2L. Notice that the product 

is known as the degrees of freedom of the interval 

2L. 

The important observation here is the exponential decay of the 
eignvalues corresponding to Indexes n larger than the inefex which is 
the degree of freedom of the interval. One concludes that the terms of 
the norm of g(x) corresponding to n > N dr are of the form of an exponent 
tlally increasing sequence, l.e., 


*n * a n C1 + e > 


bv -1 


n > N 


dr 


C7) 


a n 

One might insist that the product should converge and obtain 

A n 

a converging result. This requirement however eliminates interesting 
and useful applications. In order to counter this difficulty the 
problem will be looked at from a different point of view in the next 


-3t 


section. 


3. 


THE GREEN’S FUNCTION EXPANSION 


Observing thet ^ is the reel peat of the free apace 

Green's function one for axially synmetric fields obtains the expansion 


" Z(2n+1) J n (kx) d n (kx,) { 


f p >> > 


X-X* I ■ lx |« lx' 


<8) 


P n Cir) 

% ■ ti 

P (0) f 1 the integral equation becomes P ( tt> » (rli 


g G n J n Ckx) * F(x) , |x| < L/2 


C9i 


where 


L/2 

G n - (2n+l) J" [g(x') - (-l) n gC-x')] J n (x’) dx'» a constant &01 
0 


The problem now is that of finding some convenient numerical technique 
for a) deciding when to truncate the series to some point N and b) 
inverting the N x N matrix below 

N 

l' C n J n^ B F( V n«0» m-0,. f ,..N <H) 

n 


For small values of x i.e*. x < n the Bessel functions are small numbers 

Tt| 

and the matrix becomes ill conditioned. In the next section an estimate 
on the required number of terms will be made. 


*»*»*« 




4. THE APPROXIMATE DIMENSION OF INDEX SPACE 

In order to avrrlve at soma precise criterion of performance wc 

consider the following problem. Consider e plane wave distribution 

lit* — — 

e which results from a step function distribution X00« The dlstrl** 

but ion describes a source of angular sire a thus establishing the 

resolution limit of the reconstructed aperture. 

One obtains the one dimensional distribution in the image plane 

below 

for P - kx and u » coe 0 


F(x> 


00 — — 

- J te lktr I(k) dk x l 

—oo 


•f IX 

J e^ u A(u) du 


- 00 


For a constant source of A(u) - 1 * l u x l < cos <* 

one obtains A(u) ■ 0 , ju| > cos a 



. 2u 


sin u^ p 


u x p 


u^ - cos a 

Alternatively observing that 


F(p) 



p du 


( 12 ) 


C13) 


a4> 


and cos up- Z{2ti + 1) cos $p) j (p) P^ (u) 
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(15) 



where 


*00 - 2 O n J n (p) 


o n • 2 CO. r IP^ ( Ul ) - P n (u x )3 


a n - 2 co« ^ (2n + 1 ) (co.a) 


It is instructive to estimate the approximate dimension of the Index 
space. The distribution directly in terms of the angle 0 is 


?(9) - 1 , 0 < a 


In the space of the orthogonal functions (vh + l/2 p^tcoa©),} one has 


F(0) - 2 a Mi + 1/2 P (cos0) } 
n n 


a . P n-l* u l* I P n+l (u l^ 
n Jl v'SnHKL 


In order to find the effective dimension one compares the above with the 
flat spectrum of the delta function which Is 


6(6) - 2 


* {jtmji iP n cos a) } 


The spectrum is /Cn 




Using Parseyals theorem one has 



d 


(cos0) ^ E o n 


2 


[1 - cosa] 


( 22 ) 


Normalizing and comparing with N tarns of the delta function spectrum one 
obtains 


N~1 

S fa + 1/2) ■ I + agdl ■ tt .. .\ s s g ) . 
o V 


( 23 ) 


* 


where a Q is the amplitude of the first term 


1 - (cosa) 
~~~ 


( 24 ) 


and 1 - cosa is the quadratic content of F(8) 
One obtains 

2 2 

N ■ ■ i- - s a r 


or for large N and small a 



( 25 ) 


( 26 ) 


This result complies with our notion of resolution l,e., the 
accuracy of observation la inversely proportional with the number of 
dimensions of our observations. The following conclusion can now be 
reached. 
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1. A resolution cell of dimension 20! requires 
N - | dim.n.ion. in th« «p.rtu«. 

A consequence of this Is that If one sample at the Nyqulst rata 

X 

which for antennas Is Ax p* j corresponding to a field of •view of 
- ^ < 9 < ^ * The distension of the aperture are » 2N ^ . On the 
oth.r hand if eta »P«tur. U 1... etan the .tand.td .pertur. I < ^ 
one has to oversample In order to acquire the needed required dimension- 
ality* Oversampling however Imposed conditions on the accuracy of 
measurements and this will be examined in the next section. 

5. MEASUREMENT ACCURACY 

Let us assume that we have a flat delta function spectrum v^n + l/2 
of dimensions N which is oversampled N times a segment given by 

Y Y 

-jnr" < jx| < j , b is a constant (27) 

The bessel functions for small arguments are given by 

j„(p) 6 - , P - tat (281 

" (2n + 1) 1 1 

one then can obtain the smallest term which corresponds to the smallest 
x and largest N as follows. The amplitude of the smallest normalized 
term is 

9 2 

a o »¥ia72 i u (x)r- (1 - cosa) *'2 n+T j n C*l| (29) 

where a is obtained from eqn. 24 


- 8 - 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Wort preci*U’ly from th* Appendix and recognizing that 


1 


COB CX 


r 


» 


| a Q fiTTTTl 




l 

(2N+1) 




(30) 


In order to carry out meaningful calculations one must contain the errors 
2 

to levels below the smallest term which Is 



~ 2N 2 (2N+1) 




(31) 


where p 


kL 

2b 


For large N corresponding to small a the above equation reduces to 



2N 


kL e, 
4b N 


(32) 


The above result can also cast Into a different form by considering the 
natural degrees of freedom N q corresponding to the long aperture of length 
of L. They are 


N * 2kL 
o 


The error becomes 


w 41T 


8b N 


2n 


(33) 


( 34 ) 


Finally in tarn* of an extrpolated aperture L v daflnad by 


t* * 



one obtain* 


•t 



_i 

(2kL’) 3 




4kL' 


( 35 ) 


( 36 ) 


6. NUMERICAL EXAMPLE 

Xn ordar to lllu strata tha above extrapolation the following example 
waa considered, A target In u space of the fora 

F(u) "1 |u(* < .5 

F(u) - 0 |u| > .S 

was viewed in the image plane in an aperture 

L/2 - .810 k ~ - (2u) .810 

The aperture sustains 

N - 2.2 degrees of freedom N ■ 1.610 
o 

The target however requires at least 

N ■ /- i ■ 2 

y 1 - cosa 

or more degrees of freedom. 

Using equal sampling 11 points were observed in the image plane 
from - .810 < x < + .810. Because of the even symmetry however only 6 


points were necessary* For 12 degrees of freedom one obtains that the 

error should be below 

2 -12 
< 2.92 x 10 ^ 

Figure 2 shows the results of the extrapolation that Indicate that the 
aperture has been extended by roughly e factor of four. 

If, on the other hand, only three points were taken corresponding 
to 6 degrees of freedom the error would be 

< 1.91 x 10“ 3 

Figure 1 however shows that only a factor of 2 has been achieved. 

7. DISCUSSION OF RESULTS 

In this study one more extrapolation algorithm has been Introduced 
to the many that already exist In the literature. Questions of unique- 
ness and contamination by noise were not addressed because they are 
adequately treated In the literature. The question which was examined 
here was the following: Given that we have an aperture of length L 

having N q - 2kL degrees of freedom and an angular resolution of ^ (X 
Is the wavelength) how accurate should the measurements be In order to 
Increase the effective size to L f . In the absence of noise the results 
of our study show that the accuracy of measurements has to be at least 
as good as the term e^ f shown below 


• * 


m \\ 4 r* 


Parameter L Is an Important factor In the sampling strategy. It Ilea 
In the range 1 < b < a end It la desirable that It la as close aa It la 
possible to -L. 

This Indicates that the Important region Is the area close to the 
edges of the aperture. It la this area that has to be well surveyed In 
order to produce effective extrapolation. 

The exponential behavior shown in equation (37) Is similar to the 
one Indicated by H. Landau In the asymptotic form of his eigenvalues. 

It seems that this is an Inherent property of the extrapolation problem 
which makes It so difficult to Implement In practice. 
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APPENDIX 


BOUNDS FOR BESSEL FUNCTION 


The product of Bessel Functions is 


v (*) P+<l+2k 

J p <x) J q .(x) - Z (-1)* W V (p+g +2k 4-1) 

k! T (p+q+k+1) T(p+k+l) r(q+k+l) 

for p - q - n + 1/2 and + i (x) one has 


. 2, , tt „ , , v k y Xv 2nrfk (2nf2k-l)! 2 2n+2k+2 

Vi (') 4 Z ¥ let (2nfkn) — " , ,j2 ; 


because f 


(n+|) 


— (2n - 1) ! ! 
2 n 


j n 2 (*) 


r / - ^ ^2n+2k 

Z (-1) p 


(2n+2k+l) 1 

k! (2n+k+l>! [(2nf2k+l> ! !] 2 


- p 


2n 


{ 


I(2n+1)!!] 


+R 

[ (2n+3) 1 1 J J 


o 2 " 

(2n+l) ! J 


[1 - C 2 "* 3 ) . P 2 ] + R 
(2n+3) 
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The remainder R la a poaltlve quantity because the k+1 term is always less 
than the k term. 

The ratio of the to the ^ term is 

(2n+2k+3) p 2 K , 

\ (k+1) (2n+k+2) [(2ttf2k+3) 2 ) 


2n 


[ (2n+3) ! I) 


[1 - 


(2n+3) 


] + R 


for p < 1 
One obtains -hen 




I (,2h+l) I IJ 


a - 


2n+3 


) < 4/(x) < 


2n 


[ (2n+l) It] 


or since 

T(n + -|) 



TC2n) 


r(n) 


(2n-l) 11 
2 n 



then 

C2n+1) : ! - ^ 2n+2 ^ '• - (2n+l) <2n-l) ! ! - (2n+l) — iMl 

(n+1) ! 2 n+1 (n!) 2 1 

and since ^ 

V^T (§)" e 15 > n! >v/2^ (|)" 


A- 2 


